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ABSTRACT 



We analyze the extreme high-magnification microlensing event OGLE-2008- 
BLG-279, which peaked at a maximum magnification of A ~ 1600 on 30 May 
2008. The peak of this event exhibits both finite-source effects and terrestrial 
parallax, from which we determine the mass of the lens, Mi = O.64±O.1OM , and 
its distance, Di = 4.0±0.6kpc. We rule out Jupiter-mass planetary companions 
to the lens star for projected separations in the range 0.5-20 AU. More generally, 
we find that this event was sensitive to planets with masses as small as 0.2 M e ~ 
2 MMars with projected separations near the Einstein ring (~ 3 AU). 

Subject headings: gravitational lensing, planetary systems, planetary systems: 
formation 
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Introduction 



A complete census of planets beyond the snow line will be crucial for testing the currently 
favored core-accretion theory of planet formati on since that is t he region where this model 
predicts that giant planets form. For example, llda fc Linl (120041 ) find that gas giant planets 
around solar- type stars preferentially form in the region between the snow line at 2.7 AU and 
~ 10 AU. While radial velocity and transit searches account for most of the more than 300 
planets known to date, microlensing has the ability to probe a different region of parameter 
space that reaches far beyond the snow line and down to Earth -mass planets. Microl ensing is 
most sensitive to planets near the Einstein ring radius, which iGould fc Loebl (jl992j) showed 
lies just outside the snow line: 



Mi_ 



1/2 



AU, 



(1) 



for reasonable assumptions. This sensitivity to planets beyond the snow line is demonstrated 
by the eight published planets found b y microlensing, wh i ch range in mass from super-Earths 



to Jupiters and more massive objects (IBond et al. 



2006; Gould ct al. 2006; Bennett et al 



2008 



2004; Udals 



Dong et al. 



a et al 



200 



9; 



2005 



Gaudi et al 



Beaulieu et al. 



2008 



In high magnification microlensing events (A > 100), the images finely probe the full 
angular extent of the Einstein ri ng, making these events particularly sensitive to planets 
over a wide range of separations (IGriest fc Safizadehllf 9981 ). Additionally, because the time 
of maximum sensitivity to planets (the peak of the event) can be determined in advance, 
intensive observations can be planned resulting in improved coverage of the event, particu- 
larly given limited resources. Even when a planet is not detected, the extreme sensitivity of 
such an event can be used to put broad constraints on planetary companions. 

High magnification events are also useful because it is more likely t hat secondar y effects 
such as the finite-source effect and terrestrial parallax can be measured (jGouldll 19971 ). These 
effects can be used to break several microlensing degeneracies and allow a measurement of 
the mass of the lens and its distance. This allows us to determine a true mass of a planet 
rather than the planet/star mass ratio and a true projected separation rather than a relative 
one. Thus, in addition to being more sensitive to planets, high magnification events allow 
us to make more specific inferences about the nature of the system. 

Previous work has empirically demonstrated the sensitivity of high magnification events 
to giant planets by analyzing observed events without detected planets and explicitly com- 
puting the detection sensitivity of these events to planetary companions. Th e first high 



magn i fication event to b e analyzed in such a way was MACHO 1998-BLG-35 ( IRhie et al. 



20001 ). IRhie et al.l (120001 ) found that planets with a Jupiter- mass ratio (q = 10 3 ) were ex- 
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eluded for projected separations in units of the Einstein ring radius of d =0.37-2.70. Since 
then, many other aut hors have analyzed the planet detection sensitivity o f individual high- 
magnification events faond et alihooi baudi et allbopj IYoo et alJlipO^ : I Abe et alJliooi 
Dong et al.l 120061 ; iBatista et al.ll2009l ). In particular, prior to the work presented here, the 
most sensitive event with the broadest constrain ts on planetary c o mpanions wa s MOA 2003- 
BLG-32, which reached a magnification of 520 (lAbe et al.l 120041 ) . iDong et al.l (120061 ) found 
that this event had sensitivity to giant planets out to d < 40 

This paper presents the analysis of OGLE-2008-BLG-279, which reached a magnification 
of A ~ 1600 and was well-covered over the peak, making it extremely sensitive to planetary 
companions. In fact, as we will show, this event has the greatest sensitivity to planetary 
companions of any event yet analyzed, and we can exclude planets over a wide range of 
separations and masses. Furthermore, this event exhibited finite-source effects and terrestrial 
parallax, allowing a measurement of the mass and distance to the lens. This allows us to 
place constraints on planets in terms of their mass and projected separation in physical units. 
We begin by describing the data collection and alert process in §|2j In $3] we describe our 
fits to the light curve and the source parameters. We then go on to discuss the blended 
light and the shear contributed by a nearby star in §HJ Finally, we place limits on planetary 
companions in £j5j We conclude in §|6j 



2. Data Collection 

On 2008 May 13 (HID' = HJD - 2,450,000 = 4600.3604), the OGLE collaboration 
announced the discovery of a new microlensing event candidate OGLE-2008-BLG-279 at 
RA=17 h 58 m 36fl7 Dec=-30°22'08."4 (J2000.0). This event was independently announced by 
the MOA collaboration on 2008 May 26 as MOA-2008-BLG-225. Based on the available 
OGLE and MOA data, //FUN began observations of this event on 2008 May 27 from the 
CTIO SMARTS 1.3m in Chile, acquiring observations in both the V and I bands, and the 
next day identified it as likely to reach very high magnification two days hence. This event 
was monitored intensively over the peak by MOA, the PLANET collaboration, and many 
/iFUN observatories. Specifically, the //FUN observatories Bronberg, Hunters Hill, Farm 
Cove, and Wise obtained data over the peak of this event (see Fig. [I]). OGLE-2008-BLG- 
279 peaked on 2008 May 30 at HJD' = 4617.3481 with a magnification A ~ 1600. 



] Dong et al. ( 2006T ) also analyzed the event OGLE-2004-BLG-343, which reached a peak magnification 
of A ~ 3000. Although this is the highest magnification event analyzed for planets, sparse observational 
coverage over the peak greatly reduced its sensitivity. 
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Fig. 1. — Light curve of OGLE-2008-BLG-279 near its peak. The left panel shows the entire 
event, while the right panel shows a close-up of the peak with residuals from the point-lens 
model including finite-source effects. The black solid line shows this best-fit model. For 
clarity, the data have been binned and rescaled to the OGLE flux. 
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Because there were so many data sets, this analysis focuses on the /iFUN data from 
observatories that covered the peak of the event (/iFUN Bronberg (South Africa), Hunters 
Hill (Australia), Farm Cove (New Zealand), and Wise (Israel)) and PLANET Canopus 
(Australia) combined with the data from OGLE and MOA which cover both peak and 
baseline. We used the data from CTIO to measure the colors of the event but not in other 
analyses. Early fits of the data indicated that the Bronberg data from HJD'4617.0-4617.32 
suffer from systematic residuals that are more severe than those seen in any of the other 
data, so these data were excluded from subsequent analysis. 



The data were all reduced using difference imaging analysis (DIA; IWozniak! 120001 ) 



with the except i on of the CTIO data which were reduced using the DoPHOT package 



( ISchechter et al.l Il993l ). The uncertainties in all the data sets were normalized so that 
the x 2 /degree-of-freedom ~ 1, and we removed > 3a outliers whose deviations were not 
confirmed by near simultaneous data from other observatories. The normalization factors 
for each observatory are as follows: OGLE(1.8), MOA(l.O), /iFUN Bronberg(1.4), /iFUN 
Hunters Hill 1(2.7) and U(1.5), /iFUN Farm Cove(2.1), /iFUN Wise(3.8), PLANET Cano- 
pus(4.6), and /iFUN CTIO 1(1.4) and V(2.0). 



3. Point-Lens Analysis 

The data for OGLE-2008-BLG-279 appear to be consistent with a very-high magnifica- 
tion, A = 1570 ± 120, single lens microlensing event. We therefore begin by fitting the data 
with a point-lens model and then go on to place limits on planetary companions in §[5j In 
this section, we describe our fits to the data and address the second-order, finite-source and 
parallax effects on the light curve. 



3.1. Angular Einstein Ring Radius 

From the V- and J-band images taken with CTIO both during the peak and after the 
event, we construct a CMD of the event (Fig. [2]). We calibrate this CMD using stars 
that are also in the calibrated OGLE-III field. For the source, we measure [J, (V — I)] = 
[21.39 ± 0.09, 2.53 ± 0.01]. If we assume that the source is in the bulge and thus behind the 
same amount of dust as the clump, we can compute the dereddened color and magnitude. 
We measure the color and magnitude of the clump: [J, (V— I)] c i = [16.48, 2.71]. The absolute 
color and magnitude of the clump are [Mi, (V — I)o] c i — [—0.20, 1.05], which at a distance 
of 8.0 kpc would appear to be [I, {V - J)] 0)C i = [14.32, 1.05]. We find Aj = J d - I 0)d = 



Table I. Light Curve Fits 



Effects Fit Parameters 



JT 1111 L" 


VJL UlLdl 


T^prrocf ri p 1 
± t;l 1 Coll Idl 






+ n _ zLfil 7 34 


Uq 




(-E 


p* 


tte,e 


tte,n 




M P TP 1 1 P V 
L dl dlld-A. 


PaVQ 11 O \r 
L ObL ( X 1 1 <X -A. 


Uq 


-Av 2 


l^uays; 






l^uays; 


\VE) 






[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 




[8] 


[9] 


[10] 


[11] 


/ 








0.00 


0.00783(7) 


6.4(5) x 10" 


-4 


111.(9) 


6.6(6)xl0 4 






/ 


/ 


/ 




164.50 


0.00787(8) 


6.6(5)xl0" 


-4 


106.(9) 


6.8(6)xl0 4 


-0.15(2) 


0.02(2) 


/ 


/ 


/ 


/ 


127.97 


0.0081(1) 


-6.4(6) xlO" 


-4 


109.(9) 


6.7(6)xl0 4 


0.11(2) 


0.09(2) 


/ 


/ 






15.52 


0.00784(8) 


8.(1) xlO- 


-4 


84.(12) 


9.(1) xlO 4 


1.5(4) 


-0.3(2) 


/ 


/ 




/ 


15.51 


0.00786(8) 


-8.(1) xlO" 


-4 


84.(12) 


9.(1) xlO 4 


1.5(4) 


-0.3(2) 


/ 




/ 




166.40 


0.00787(8) 


6.9(6)xl0" 


-4 


101.(8) 


7.2(6) xlO 4 


-0.16(2) 


0.03(2) 


/ 




/ 


/ 


129.59 


0.0081(1) 


-6.9(5)xl0" 


-4 


102.(8) 


7.1(6)xl0 4 


0.11(2) 


0.11(3) [ 



Note. — The first 4 columns indicate which effects were included in the point-lens fit. The Ax 2 improvement for each fit 
(col. 5) is given relative to the best-fit including finite-source effects but without parallax. There are 5731 data points in the 
fit light curve. The numbers in parentheses indicate the uncertainty in the final digit or digits of the fit parameters. 
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Fig. 2. — Calibrated Color-Magnitude Diagram (CMD) constructed from the CTIO and 
OGLE data. The square indicates the centroid of the red clump, the open circle shows the 
blended light, and the solid circle indicates the source. The small black points are field stars. 
The error bars are shown but are smaller than the size of the points. 
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16.48 - 14.32 = 2.16 and E(V - I) = (V - I) d - (V - 7) 0jC i = 1-66. We then calculate the 
dereddened color and magnitude of the source to be [I, (V — i)] = [19.23, 0.87]. 

The angular Einstein ring radius can be determined by combining information from the 
light curve and the color-magnitude diagram (CMD). Finite source effects in the light curve 
enable us to determine the ratio of the source size, 9*, to the Einstein radius, 9 E : 

P* = 9J9 E . (2) 

We can then estimate 9+ from the color and magnitude of the source measured from the 
CMD, and solve for 9 E . 



3.1.1. Finite-Source Effects 



If the source passes very close to the lens star, finite-source effects will smooth out the 
peak of the light curve and allow a measurement of the source size p*. Although finite- 
source effects are not obvious from a visual inspection of the light curve, including them 
yields a dramatic improvement in \ 2 . In order to fit for finite source effects, we first estimate 
the limb-darkening of the source from its color and magn itude. We combine th e color and 
magnitude of the source with the Yale-Yonsei isochrones (jDemarque et all 12004 ). assuming 
a distance of D s = 8 kpc and solar metallicity, to estimate T eS = 525 0K an d log q = 4.5. We 
use these values to calculate the limb-darkening coefficients, u, from I Clare tl (120001 ). assuming 
a microturblent velocity of 2 km/s. We calculate the linear limb-darkening parameters IV 
and Tj using T = 2u/(3 — u) to find T v = 0.65 and Tj = 0.47. We use these values in our 
finite-source fits to the data. We find that a point-lens fit including finite-source effects is 
preferred by A% 2 of 2647.85 over a fit assuming a point source. We search a grid of uq and 
p* near the minimum t o confirm that thi s is a well constrained result. We use zq = ito/p+ as 
a proxy for p* following lYoo et al.l (120041 ) . The resultant x 2 niap in the Uq-z plane is shown 
in Figure [3l Our best-fit value for p* is 6.6 ± 0.6 x 10~ 4 . For this value of p*, z is almost 
unity, indicating that the source just barely grazed the lens star. The other parameters for 
our best-fit including finite-source effects are given in Table [H 



3.1.2. Source Size 



We convert the dereddened color and magnitude of the source to (V— if) using 



Bcsscl 



k Brett 



(119881 ). and combine them with the surface brightness relations in iKervella et al.l (120041 ) to 
derive a source size of 0* = 0.54 ±0.4 pas. The uncertainty in 9* comes from two sources: the 
uncertainty in the flux and the uncertainty in the conversion from the observed (V — I) color 
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Fig. 3. — x 2 contours as a function of impact parameter, u , and z = u / p* where p* = 0+/0e 
is the normalized source size. The best fit is marked with a plus sign. 
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to surface brightness. The uncertainty in the flux (i.e. the model fit parameter f s j) is 8.5%, 
and we adopt 7% as the uncertainty due to the surface brightness conversion. From equation 
(J2J), we find that #e = #*/p* = 0.81 ± 0.07 mas. We also calculate the (geocentric) proper 
motion of the source: /z geo = # E /t E = 2.7 ± 0.2 mas/yr. Because the peak flux (oc f s ,i/p±) 
and source crossing time (p*£e) are both essentially direct observables, and so are well con- 
strained by the light curve, the fractional uncertainty in Be and /i geo are comparable to the 
fractional uncertainty in 9*. This result is generally applicable to point-lens/finite-source 
events and is discussed in detail in the Appendix. 



3.2. Parallax 



Given that we have a measurement for #e, if we can also measure microlens parallax, 
7Te, we can combine these quantities to derive the mass of the lens and its distance. The 
mass of the lens is given by 



Mi 



K = 



Its distance D\ is 



kit e 
1 AU 



AG 
c 2 AU 



8.14 



mas 



M P 



D. 



7TZ = 7T S 



TTreb 



(3) 



(4) 



where 717 is the parallax of the lens, n s = 0.125 mas is the parallax of the source (assuming 
a distance of D s = 8 kpc), and 7r rel = 9eti e . 

Microlens parallax is the combination of two observable parallax effects in a microlensing 
event. Terrestrial parallax occurs because observatories located on different parts of the 
Earth have slightly different lines of sight toward the event and so observe slight differences 
in the peak magnific ation and in the timing of the peak, descr ibed by the parameters uq 
and to, respectively (IHardy fc Walker! Il995l ; iHolz fc Waldl 119961 ). Orbital parallax occurs 
because the Earth mo ves in its orbit during the event, again, changing the apparent line 
of sight. iGouldl (119971 ) argued that one might expect to measure both finite-source effects 
and terrestrial parallax in extreme high-magnification events. We fit the light curve for both 
of the sources of parallax, including finite-source effects. Fitting for both kinds of parallax 
simultaneously yields a A% 2 improvement of 165 (see Tabled]). We find 7Te = (vte,e, tte.n) = 
(—0.15 ± 0.02,0.02 ± 0.02), where 7Te,e and 7Te,n are the projections of 7Te in the East and 
North directions, respectively. 



Smith et al.l (120031 ) showed that for orbital parallax and a constant acceleration, uq has 



a si gn degenerac y. This degeneracy may be broken if terrestrial parallax is observed (see 
also IGouldl 12004 ) . In the fits described above, we assumed u > 0. We repeat the parallax 
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fit fixing uq < 0. We find that the +uq solution is preferred over the — uq case by A% 2 = 37 
(see Tabled]). 

We perform a series of fits in order to isolate the source of the parallax signal, i.e. 
whether it is primarily due to orbital parallax or terrestrial parallax. We first fit the light 
curve for orbital parallax alone and then fit for terrestrial parallax alone. The results are 
given in Table [TJ For +uo, the orbital parallax fit gives (vte,e, tte,n) = (1-5 ± 0.4, —0.3 ± 0.2) 
and a Ax 2 improvement of ~ 16 over the finite-source fit without parallax. In contrast, 
the +v,o fit for terrestrial parallax alone yields A% 2 = 166 and (7Te,E) 7Te,n) = (—0.16 ± 
0.02,0.03 ± 0.02). While the orbital and t errestrial parallaxes ar e nominally inconsistent 



at more than 3a, from previous experience (IPoindexter et al.ll2005l ) we know that low- level 



orbital parallax can be caused by small systematic errors or xallarap (the orbital motion of 
the source due to a companion), so we ignore this discrepancy. From the A% 2 values, it is 
clear that terrestrial parallax dominates the microlens parallax signal in this event, so any 
spurious orbital parallax signal does not affect our final results. 

We also confirm that the terrestrial parallax signal is seen in multiple observatories, and 
thus cannot be attributed to systematics in a single data set. To test this, we repeat the fits 
for parallax excluding the data from an individual observatory. If a data set is removed and 
the parallax becomes consistent with zero, then that observatory contributed significantly to 
the detection of the signal. Using this process of elimination, we find that the signal comes 
primarily from the MOA and Bronberg data sets. 

Given the results of these various fits, we conclude that the best fit to the data is for the 
+uq solution, and we include both forms of parallax for internal consistency. Combining this 
parallax measurement with our measurement of 8e from §3.11 we find Mi = 0.64 ± 0.1M Q 
and Di = 4.0 ± 0.6 kpc (7r re i = 0.13 ± 0.02 mas) using equations ([3]) and (J1J. 



4. The Blended Light 

The centroid of the light at baseline when the source is faint is different from the centroid 
at peak magnification, indicating that light from a third star is blended into the PSF. The 
measured color and magnitude of blended light are [I, (V — /)]& = [17.21 ± 0.01, 2.32 ± 0.02]. 
Stars of this magnitude are relatively rare, and so the most plausible initial guess is that the 
third star is either a companion to the source or a companion to the lens. If the former, we 
can use the values of Aj and E(V — I) we found above to derive the intrinsic color of the 
blend: [I, (V — I)]o,b = [15.05, 0.66]. This assumes that the blend is in the bulge at a distance 
of 8 kpc, giving an absolute magnitude of Mj^ = 0.53 and M Vtb = 1.19. Figure @] shows 
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V-l 



Fig. 4. — Possible absolute magn itudes and colors for the blend plotted with Yale-Yonsei 
isochrones (IDemarque et al.ll2004l ). The isochrones plotted are the Y 2 isochrones for solar 
(thick) and sub-solar metallicities (thin) for populations 1 (dotted), 5 (dot-dashed), and 10 
Gyr old (solid). The dashed line shows the color and magnitude of the blend for a continuous 
distribution of distances assuming a dust model that decreases exponentially with scale 
height. The square shows the absolute magnitude and color of the blend assuming it has the 
same distance (8 kpc) and reddening as the clump. The plus sign, diamond, and triangle 
show the absolute magnitude and color using our simple dust model and distances of 2, 4, 
and 6 kpc, respectively. If the blend is a companion to the lens, it would be at a distance of 
4 kpc (diamond). 
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this point (open square) compared to solar (Z=0.02) and sub-sol ar metallicity (Z=0.001) 
Yale-Yonsei isochrones at 1, 5, and 10 Gyrs (IDemarque et al.l 12004 ). These isochrones show 
that the values of [My, (V — I)o]b may be consistent with a sub-giant that is a couple Gyr 
old, but a more precise determination of age is not possible since the age is degenerate with 
the unknown metallicity of the blend. 

If the blend is a companion to the lens, however, it lies in front of some fraction of the 
dust. In order to derive a dereddened color and absolute magnitude to this star, we need a 
model for the dust. We explore this scenario using a simple model for the extinction that is 
constant in the plane of the disk and decreases exponentially out of the plane with a scale 
height of H = 100 pc: 

, -Dsinb^ 

1 — exp 



Md) = K x 



Hn 



(5) 



where D is the distance to a given point along the line of sight, h is the Galactic latitude, 
and K\ is a constant. We can solve for K\ by substituting in the value of Aj that we find 
for the source at 8 kpc. We then model the selective extinction in a similar manner: 



E(V -I) = K 2 



1 — exp 



— D sin b 



(6) 



and solve for K 2 using the value of E(V-I) calculated for the source at 8 kpc. From equations 
|5]and[6l we can recover the intrinsic color and magnitude of the blend assuming it is at various 
distances. In Figure HJ we plot a point assuming the blend is at the distance of the lens, 
4.0 kpc. By interpolating the isochrones and assuming a solar metallicity, we find that the 
blend is consistent with being a 1.4 M sub-giant companion to the lens with an age of 3.8 
Gyr. For comparison, we also plot a line showing how the inferred color and magnitude of 
the blend vary with the assumed distance. 



4.1. Astrometric Offset 

From the measured blend flux, one can determine the astrometric offset of the source 
and blend by comparing the centroid of light during and after the event. At a given epoch, 
the centroid is determined by the ratio of the flux of the blend to the sum of the fluxes of 
the source and lens. That ratio depends on the magnification of the source. Thus, if we 
know the magnification of the source at two different epochs and the intrinsic magnitude 
of the source and the blend, we can solve for the separation of the lens and the blend. We 
find A6 = 153 ± 18 mas. Given this offset, we will show below that based on the lack of 
shear observed in the light curve, the blended light cannot lie far in the foreground and thus 
cannot be the sub-giant companion to the lens hypothesized above. 
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4.2. Search for Shear 

Because all stars have gravity, if the blend described above lies between the observer and 
the source, it will induce a shear 7 in the light curve. We can estimate the size of the shear 
using the observed astrometric offset and assuming that the blend is a 1.4M companion to 
the lens. 

7 AO 2 ~ AO 2 ' 

= 6.2 x ,0- (-?=!-) (*) (-*»-) - 2 . (7, 
V 0.13 mas/ \IAM & ) \ 153 mas/ v ; 

Using the lex upper limit on the separation (171 mas), we find a minimum shear of 7 = 
4.9 x 10~ 5 if the blend is a companion to the lens. To determine if this value is consistent 
with the light curve, we perform a series of fits to the data using binary-lens models that 
cover a wide range of potential shears. The effect of the shear is to introduce two small 
bumps into the light curve as the small binary caustic crosses the limb of the source, and 
this is indeed what we see in the binary-lens models we calculate. 

Because the separation between the lens and a companion is large (B = AO/ 9-^ ^> 1), 
the shear can be approximated as 7 ~ Q/B 2 , where Q = M^/Mi is the mass ratio of the 
companion and the lens. This reduces the number of parameters that need to be considered 
from three to two: 7 and a, the angular position of the blend with respect to the motion 
of the source. We use a grid search of 7 and a to place limits on the shear. For each 
combination of 7 and a, we generate a binary light curve in the limit B ^> 1 that satisfies 
Q = 7-B 2 and fit it to the data using a Markov Chain Monte Carlo with 1000 links. We bin 
the data over the peak to reduce computing time. We compute the difference in \ 2 between 
the binary model and the best-fit finite-source point-lens model. Figure shows the results 
of the grid search over-plotted with the upper and lower limits on the shear assuming the 
blend is a companion to the lens. From this figure, we infer that a shear of 6.2 x 10~ 5 is 
inconsistent with our data since it is in a region where the fit is worse by A\ 2 > 36. 

The two minima in the x 2 map at 7 ~ 10~ 4 ,a = tt/2,tt are well-defined but appear 
to be due to a single, deviant data point. Fits to the data with these binary models show 
improvement in the fit to this data point, but the residuals from these fits for the other 
data points are large and show increased structure. Thus, we believe these minima to be 
spurious and conclude that the maximum shear that is consistent with our data (A% 2 < 9) 
IS 7max = 1.6 x 10~ 5 . 

Since we have ruled out the scenario where the blend is a companion to the lens, we 
need to ask what possible explanations for the blend are consistent both with 7 max and with 
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Fig. 5. — Shear as a function of a (angular position with respect to the motion of the source). 
Open symbols indicate an improved x 2 compared to the finite-source point-lens fit. Filled 
symbols indicate a worse fit. The magnitude of Ax 2 is indicated by the color legend shown. 
The solid line indicates our calculated value for the shear assuming the blend is at the same 
distance as the lens. The shaded area shows the 1 a limits on this value from the uncertainty 
in the centroid of the PSF (see text). 
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the observed color and magnitude. Given 7 max , we can place constraints on the distance to 
the blend, D b , for a given mass. The distance is given by 

D b = -, (8) 

7(A#) 2 

where n b = tt s + n Tcl b = n s H — — (9) 

KMf, 

If we assume M b = 1M , 7 = 7 max , and use previously stated values for the other parameters, 
we find D b > 5.8 kpc. A metal-poor sub-giant with this mass located at or beyond this 
distance would be consistent with the observed color and magnitude of the blend given the 
simple extinction model described above. However, other explanations are also possible. For 
example, if the mass of the blend were decreased, n b would increase, and a slightly closer 
distance would be permitted. Thus, we cannot definitively identify the source of the blended 
light. However, given that 7 max is very small, we can ignore any potential shear contribution 
in later analysis. 



Limits on Planets 



We use the method described by IRhie et al.1 (120001 ) to quantify the sensitivity of this 
event to planets. This approach is used for events such as this one for which the residuals are 
consistent with a point-lens. Rather than fitt i ng bi nary models for planetary companions 
to our data as advocated by lGaudi fc Sackettl (120001 ) . we generate a binary model from the 
data and fit it with a point-lens model. When the single-lens parameters are well constrained 
(as is the case with O GLE-2008-BLG-27 9) , t hese two approac hes are essentially equivalent 
(see the discussion in iGaudi et alJl2002l and iDong et al.1 120061 ). We create a magnification 
map assuming an impact parameter, d, and star/planet mass ratio, q, using a lens with the 
characteristics from our finite-source fit. The method for creatin g the magnification map 



is described in detail in IDong et al.1 (120061 ) and IDong et al.l (120091 ) . For each epoch of our 
data, we generate a magnification due to the binary lens assuming some position angle, a, 
of the source's trajectory relative to the axis of the binary and assign it the uncertainty of 
the datum at that epoch. As in §4.2[ we use binned data for this analysis. 

For q = 10 -3 , 10~ 4 , 10~ 5 , and 10~ 6 we search a grid of d, a and c ompute the Ay 2 . 
Base d on the systematics in our data, we choose a threshold Ax^ in = 160 (IGaudi fc Sackett 



20001 ). For Ax 2 > Ax 2 ^ , the fit is excluded by our data, and we are sensitive to a planet 



of mass ratio q at that location. We repeat the analysis using unbinned data for a small 
subset of points and confirm that the A% 2 for fits with the unbinned data is comparable 
to fits with binned data. Figure [6] shows the sensitivity maps for four values of q. These 
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-1 1 -1.0 -0.5 0.0 0.5 1.0 

d cos a d cos a 



Fig. 6. — Planet sensitivity as a function of distance from the lens in units of Einstein 
radii. The white/black circle indicates the Einstein ring (d = 1). The mass ratios and 
corresponding planet masses are indicated on each plot. The colors indicate the A% 2 that 
would be caused by a planet at that location. 
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Fig. 7. — Detection efficiency map in the (d, q) plane, i.e. projected separation in units 
of #e and planet-star mass ratio. The contours show detection efficiencies of 0.99, 0.90, 
0.75, 0.50, 0.25, and 0.10 from inside to outside. The inner spike is due to resonant caustic 
effects at the Einstein ring. The upper and right axes translate (d, q) into physical units 
(r±,m p ), i.e. physical projected separation and planet mass. The vertical solid line shows 
the position of the snow line for this star. The dotted line shows Kepler's sensitivity to 
planets around the lens star assuming my = 12. The cutoff in separation (d ~ 0.6) occurs 
where a planet's orbital period is equal to Kepler's mission lifetime of 3.5 yrs. The dashed 
line shows the sensitivity limit for radial velocity observations with 1 m/s precision. The dot- 
dashed line shows the sensitivity limit for a space-based astrometry mission with precision 
of 3 /zas assuming the star is at 10 pc. 
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maps show good sensitivity to planets with mass ratios q = 10~ 3 , 10~ 4 , and 10~ 5 and some 
sensitivity to planets with q = 10 -6 . For our measured value of Afj = O.64M , a mass ratio 
of q = 10~ 3 corresponds to a planet mass m p = 0.67Mj up and a mass ratio of q = 10~ 6 
corresponds to m p ~ 2MMars- The results bear a striking resemblance to the hypothetical 
planet sensitiv ity of the Ama* ~ 3000 event OGLE-2004-BLG-343 if it had been observed 



over the peak (IDong et al.ll2006l ). In particular, this event shows nearly uniform sensitivity 
to planets at all angles a for large mass ratios. The hexagonal shape of the sensitivity map 
is the imprint of the difference between the magnification maps of pl anetary-lens mode ls and 
their corresponding single- lens models (see upper panel of Fig. 3 in IDong et al.l (120091 ) ). 



Figure [7] shows a map of the planet detection efficiency for this event. The efficiency is 
the percentage of trajec tories, a, at a given mass ratio and separation that have A% 2 > Ax^ in 
( iGaudi fc Sackettl l2000h . The efficiency contours are all quite close together because of the 
angular symmetry described above for the planet sensitivity maps. Because we measure the 
distance to the lens, we know the projected separation, rj_, in physical units: 



r ± = d0 E Di. 



(10) 



Since we know Mi, we also know the planet mass, m p = qMi. We can rule out Neptune- mass 
planets with projected separations of 1.5-7.2 AU (d = 0.5-2.2) and Jupiter-mass planets 
with separations of 0.54-19.5 AU (d = 0.2-6.0). We are also able to detect Earth-mass 
planets near the Einstein ring, although the efficiency is low. The region where this event is 
sensitive to giant planets probes well beyond the s now line of this star, which we estimate 
to be at 1.1 AU assuming a snow = 2.7AU(M*/M ) 2 (llda fc Linll2004l ). The observed absence 
of planets, especially Neptunes, immediately beyond the snow line of this star is interesting 
given that core-accretion th eory predicts that N e ptune-mass plan ets should preferentially 
form around low-mass stars (ILaughlin et al.ll2004l ; llda fc Lid 120051 ). 



It is also interesting to consider how the sensitivity of this event to planets compares to 
the sensitivity of other planet-search techniques. Obviously, because of the long timescales 
involved, most transit searches barely probe the region of sensitivity for this event. As a 
space-based mission, the Kepler satellite has the best opportunity to probe some of th e 
microlensing parameter space using transits. Using equation 21 from lGaudi fc Winn! (120071 ). 
we can estimate Kepler's sensitivity to transits around this star: 



m r 



0.22 ( S -B) 3/2 



10 



1 AU 



X 3/4 

J lO a3(m ^- 12) Mi 



Earth j 



where {S/N) is the signal-to-noise ratio, a is the semi-major axis of the planet, and my is 
the apparent magnitude of the star. We have assumed that the density of the planet i s the 
same as the density of the Earth and the stellar mass-radius relation = kM®' 8 (ICox 
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2000l ). Kepler is also limited by its mission lifetime of 3.5 yrs. Fo r periods longer tha n this, 



it becomes increasingly unlikely that Kepler will observe a transit (lYee fc Gaudill2008l ). This 
limits the sensitivity to planets within ~ 2 AU where the period is less than the mission 
lifetime. These boundaries are plotted in Figure [7J 

For comparison, we can also estimate the sensitivity of the radial velocity technique to 
planets around a star of this mass assuming circular orbits and an edge-on system. Radial 
velocity is sensitive to planets of mass 

where ctrv is the precision, and A" is the number of observations. The limit of radial velocity 
sensitivity is plotted in Figure [7J as a function of separation assuming a precision of 1 m/s. 
Additionally, we can consider how this microlensing event compares to the sensitivity of a 
space-based astrometry mission with microarcsecond precision (a a = 3 //as): 

m „ = 6 .4 (j^.) ( S JH) (JL) - 1 ' 2 ( « V' (_L.) M E „ t , (13) 

We assume circular face-on orbits. We show the limiting mass as a function of semi-major 
axis in Figure [7J for 3 //as precision. While these contours encompass a large region of the 
parameter space, they do not take into account the time it takes to make the observations, 
which increases with increasing semi-major axis. Furthermore, we only expect this kind of 
astrometric precision from a future space mission, whereas this event shows that microlensing 
is currently capable of finding these planets from the ground. This discussion shows that 
microlensing is sensitive to planets in regions not probed by transits and radial velocity and 
will be particularly important for finding planets at wide separations where the periods are 
long. For example, for semi-major axis a = 4AU (near the maximum sensitivity shown in 
Fig. [7J, the period is P ~ 10 yr. 



6. Summary 



The extreme magnification microlensing event OGLE-2008-BLG-279 allowed us to place 
broad constraints on planets around the lens star. Even with a more conservative detection 
threshold (Ax 2 > 160), this event is more sensi tive than any pr eviously analyzed event 
(the prior record holder was MOA-2003-BLG-32; lAbe et al.l 12004 ) . Furthermore, because 
we observe both parallax and finite-source effects in this event, we are able to measure the 
mass and distance of an isolated star (M; = 0.64 ± 0.10M G , Di = 4.0 ± 0.6kpc). Using these 
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properties of the lens star, we convert the mass ratio and projected separation to physical 
units. We can exclude giant planets around the lens star in the entire region where they are 
expected to form, out beyond the snow line. For example, Jupiter-mass planets are excluded 
from 0.54-19.5 AU. Events like this that can detect or exclude a broad range of planetary 
systems out beyond the snow line will be important for determining the planet frequency at 
large separations and constraining models of planet formation and migration. 

We acknowledge the following support: NSF AST-0757888 (AG,SD,JCY); NASA NNG04GL51G 
(DD,AG,RP); Polish MNiSW N20303032/4275 (AU); Korea Astronomy and Space Science 
Institute (B-GP,C-UL);Creative Research Initiative Program (2009-008561) of Korea Science 
and Engineering Foundation (CH). 



7. Appendix: Uncertainty in 9* /i, and 9 E 

In the present case, the fractional errors in 9*, fi, and 9 E are all very nearly the same, 
although for somewhat different reasons. Since the same convergence of errors is likely to 
occur in many point-lens/finite-source events, we briefly summarize why this is the case. We 
first write (generally), 

0, = yfZ/Z 

where f s is the source flux as determined from the model, and Z is the remaining set of 
factors, which generally include the surface brightness of the source, uncertainties due to the 
calibration of the source flux, and numerical constants. Next, we write 

— He — — — a — Hi - \ t 

where /grand = fs/ P and t* = pt E . We note that for point-lens events with strongly detected 
finite source effects, t* and / gra nd are quasi-observables, and so have extremely small errors. 
For example, if Uq = 0, then 2t* is just the observed source crossing time while 2/ gran d[l + 
(37r/8 — 1)T] is the observed peak flux. Even for uq ^ 0, these quantities are very strongly 
constrained, with errors o"/ grand = 0.4% and a u = 0.3% in the present case. Since the 
errors in f s and Z are independent, the fractional errors in 9*, p, and 9e are each equal to 
[(1/ '4) iofj f s ) 2 + (a z /Z) 2 } 1 / 2 . In the present case, GfJ f s is given by the fitting code to be 
8.5%, while we estimate cr z /Z to be 7%, and therefore find a net error in all three quantities 
(0„0e, and fi) of 8%. 
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